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Abstract 



We examine correlations of energy density induced by initial state fluctuations, which are lo- 
calized in both transverse and longitudinal extent. The hotspots are evolved according to hy- 
drodynamics in a background which includes radial flow. Two-point energy density correlations 
from these hotspots are computed as a function of the difference in azimuthal angle and rapidity. 
Such localized perturbations occur naturally in the theory of hydrodynamic fluctuations and may 
provide insight into some features of the two-particle correlation data from RHIC and the LHC. 



1. Introduction 

Two-particle correlation measurements in heavy ion collisions (U O [3] [4) contain many quali- 
tatively different features than those in pp collisions; the hydrodynamic nature of the quark-gluon 
plasma has a significant effect on correlations. We are interested in the impact of local fluctu- 
ations on two-particle correlation data from RHIC and LHC. In contrast to typical initial state 
fluctuations, which are assumed to be nearly boost invariant, here the fluctuations are localized 
in both transverse and longitudinal extent. In [5 ] it was shown that initial state hotspots placed 
at a particular position in the transverse plane could reproduce many of the features of the two- 
particle azimuthal correlations. Our approach follows this work closely, but allows us to also 
examine the rapidity extent of the correlations. 

Point-like fluctuations occur naturally due to thermal noise in hydrodynamics. In (6) it was 
shown that local hydrodynamic fluctuations occurring throughout the evolution of the system 
could induce long range rapidity correlations. These proceedings are a first step toward the 
extension of these results to include the angular dependence of the correlation functions. 

2. Background flow 

In the same manner as (H, we consider an exact solution to the equations of hydrodynamics 
which incorporates radial flow. The solution we use is the Gubser flow (THS). The equations of 
hydrodynamics are solved exactly for a fluid which is at rest in the coordinate system 



The angle (p is the usual azimuthal angle around the beam line, and £ = tanh" (z/t) is the space- 
time rapidity. The "conformal coordinates" p, 6 are related to the usual proper time r = V? 2 - z 2 , 
and radial distance from the beam line r = ^x 2 + y 2 by the coordinate transformations 



ds 2 = r(p, 6) [-dp 2 + cosh 2 p (d6 2 + sin 2 6dcp 2 ) + df] . 



(1) 
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1 - (qr) 2 + {qrf 
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These coordinate transformations contain a free parameter q which has dimensions of energy; it 
allows one to specify a scale which is related to the transverse size of the system R ~ q~ l . In the 
limit q — > 0, the system has infinite transverse extent and the Bjorken solution is recovered. The 
radial velocity of the fluid is given by 



r l + (qr) 2 + (qr) 2 ' KJ 

The Gubser flow is an idealization of the quark-gluon plasma formed in a heavy ion collision. 
In particular, the solution requires azimuthal symmetry around the beam line and conformal 
symmetry. The former is restricts us to central collisions only; the latter is only approximately 
true for T » T c , where T c is the crossover temperature. 



3. Perturbations of the flow 

We add small perturbations of energy density and fluid velocity to the background flow, 

s(p, 6, <f>,& = ^ [1 + 6(p, 0, 0, 0] , (4) 

1/^,0,0,0= l - (l,n f (p,M,f)). (5) 

The background energy density £o is found by solving the equations of hydrodynamics in the 
absence of the perturbations. 

To reduce the linearized hydro equations to a set of ordinary differential equations, we employ 
a spherical harmonic expansion in the coordinates, and a Fourier transform in the longitudinal 
direction^ 

— Ys'ipwae^yv, (6) 
u (p,e,<p,o = d g f -^ Yu l s (p)Y lm (e,<p)e ik t, (i) 
uHp,0,<p,O = -^Y-d* f ^-Yd^Yue,^, (8) 

sin z 6 J (2tt) 
J (2*0 ^ ^ 

There are then three coupled equations for the perturbations 5 l , u l s ,v^, which are given in (8). 
When viscosity is neglected, the equations can be solved analytically. We assume an initial 
condition at conformal time p = po such that the energy density is initially localized, 

( 2 + 6l- 200 cos(0 - O ) + f ) 
6(po,0,<t>,& = Aexp °- 2£2 - -L (10) 

and all fluid velocities vanish. The solutions are then transformed back to the more familiar r, r 
coordinates. We plot the time evolution of an example initial state hotspot in Fig. [T] 



^ote that we have reduced u e and to the gradient of a single scalar function. In fact, there is one additional 
divergenceless vector component to the velocity, but it does not couple to the equations for the energy density. In this 
work, we will present results for the correlations of energy density, so we do not need to consider the vector piece. 
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Figure 1: Plots of the evolution of a hotspot. The columns show the profile of energy density at three different times. The 
upper row shows the evolution in the transverse plane (after an integration over rapidity), while the lower row shows the 
evolution in the r - % plane (after an integration over azimuthal angle). 

4. Correlations 

The solution to the differential equation for the over-density 6 at a final time Tf is a function 
of three coordinates r, cp, Usually, to translate this coordinate space result into momentum 
space (relevant for comparison with experiment), a Cooper-Fry e freeze-out is performed which 
involves integration over all coordinates. We will defer a full freeze out analysis to future work 
and here simply integrate only over the radial distance at constant Tf. Such an integration would 
arise from the Cooper-Fry e prescription assuming that matter at different radial positions freezes 
out at a fixed proper time Tf. As such, we define 

S((f>^)= rdrS(T f9 r,<f> 9 &. (11) 
Jo 

We compute the correlation function by multiplying two such solutions and integrating over 
the shift of azimuthal angle and rapidity. In this way, we are integrating over all hotspots placed 
at varying initial angles and rapidities. As in [5 ], for simplicity, the initial radial position of the 
hotspot, ro, is fixed and not integrated. The correlator, 

r*2n poo 

C ss {T f , A<p, A£> = dx\ dr, 6(r f , ft - X &- n)S(r f , ft - - rj). (12) 

JO J-oo 

is thus a function of the angular difference A<p = <p\ - cp2 and rapidity difference A£ = %\ - 
The results are shown in Fig. [2] 
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Figure 2: Plot of the correlation function of energy density as a function of the (coordinate space) angular difference 
and spacetime rapidity difference. To make this plot, we summed harmonics up to / = 25. The parameters used were 
to = 0.5 fm/c, Tf = lOfm/c, q~ l = 4.3 fm, I = 0.1, A = 1/2. Following | 5 1, the initial radial position of the hotspot was 
r = 4.1 fm. 



5. Summary 

We have examined the energy correlations induced by an initial hotspot localized in both 
transverse and longitudinal extent. Because the hotspot propagates according to hydrodynamics, 
a sound front is evident in Fig. [2] In addition to the sound front, we notice a peak at zero sepa- 
ration which falls off more slowly in the rapidity direction than in the angular direction. These 
properties may have some implications for near side peak observed in two-particle correlation 
measurements at RHIC and LHC. 

Rather than introducing the initial source by hand, this approach should be extended to in- 
clude stochastic "noise" in order to make this approach consistent with the theory of hydrody- 
namic fluctuations (6). It is also necessary to perform a proper freeze-out to translate the cor- 
relations in coordinate space to those in momentum space (which are observed in experiment). 
These investigations are currently underway. 
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